Bulk and Surface Nanoscale Hole Density Inhomogeneity in HgBa 2 Cu04 + 5 and 

Bi 2 Sr2CaCu20 8+( 5 Cuprates 
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It is well established that the hole density in the prototypical superconductor La2- a; Sr :E Cu04 is 
very inhomogeneous due to Sr-dopant induced disorder. On the other hand, the hole distribution 
in HgBa2Cu04+{ and Bi2Sr2CaCu2 0g+{ doped by interstitial oxygen is believed to be much more 
uniform. Recent nuclear magnetic resonance measurements indicate, however, that the charge in- 
homogeneity in HgBa2CuC>4+<s is close to that in La2-xSra;Cu04. Calculations performed in the 
present paper confirm this observation. We also show that the charge inhomogeneity is most pro- 
nounced at the surface layer that can be probed by scanning tunneling microscope. Our simulations 
demonstrate that, despite having similar amplitudes of charge inhomogeneity, the hole mean free 
path in HgBa2Cu04+i is substantially longer than that in La2- a; Sr :! ;Cu04. The screening of the 
Coulomb repulsion in HgBa2Cu04+«5 is also stronger. These two reasons might explain the difference 
in the superconducting critical temperatures between these two compounds. 

PACS numbers: 74.62.En, 74.72. Gh, 76.60.Gv 



I. INTRODUCTION 

The nanoscale electronic disorder is a long standing 
problem in physics of cuprates. This problem has many 
aspects, among which the most important one is the 
influence of disorder on critical temperature T c . Since 
the energy scale associated with pairing mechanism in 
CuC>2 plane is believed to be universal, different out-of- 
planc defects arc expected to influence T c differently. 1 ^ 
Another important aspect is the nanoscale inhomogene- 
ity of the local density of states (DOS) measured by 
scanning tunneling microscope (STM). This effect was 
mostly studied 3 - - — in Bi 2 Sr 2 CaCu208+5. In this paper, 
we address theoretically the problem of the nanoscale 
charge inhomogeneity in cuprates, and a related issue of 
the hole mean free path. We focus on two single layer 
families of high-T c superconductors, La2- 2; Sr : ;;Cu04 and 
HgBa2Cu04+a, in which the density inhomogeneity has 
been quantified by nuclear quadrupole resonance (NQR) 
experiments^^ We also calculate the charge inhomogene- 
ity in the surface layer of Bi 2 Sr2CaCu208+5- 

The NQR measures the energy splitting of nuclear lev- 
els induced by an electric field gradient at the nucleus. 
A hole in the 3d-shell of the Cu ion gives a dominant 
contribution to the field gradient at the Cu nucleus. The 
Cu NQR frequency in cuprates is therefore very sensi- 
tive to the doping level and is directly proportional to 
the local hole concentration. The contribution of the 
3d-holc to the field gradient is significantly compensated 
by contributions of holes located at nearby oxygens^ 
therefore the slope of the NQR frequency versus doping 
varies among different cuprate families. 

It is widely accepted that the hole density distribution 
within the Cu02 layer in La2- a; Sr a ;Cu04 is very nonuni- 
form, as it has been clearly demonstrated by measure- 
ments of 63 Cu NQR spectra^. The observed broad NQR 
spectrum^ unambiguously indicates a very inhomoge- 



neous hole density profile in the bulk of the sample. The 
inhomogeneity is due to the doping mechanism, where Sr 
substitutions of La ions create an effective Coulomb de- 
fect very close to the conducting Cu02 plane. In a recent 
paper,— we have performed the Hartrcc-Fock simula- 
tion of the charge density distribution in La2- 2: Sr 2: Cu04, 
which shows a very inhomogeneous charge density pro- 
file at the nanometer scale, and reproduces the observed 
NQR lincshapes. 

On the other hand, it is generally believed that hole 
density distribution in compounds doped by interstitial 
oxygens, such as HgBa2Cu04+a and Bi2Si'2CaCu208+5, 
is much more uniform than in La2- a: Sr a ;Cu04. This is 
because the distance from the interstitial oxygen to the 
C11O2 layer is typically larger than the Sr-layer distance 
in La2- a; Sr a ;Cu04. However, Cu NQR measurements in 
the single layer HgBa2Cu04+#£'ii as well as in the bi- 
layer HgBa2CaCu206+^ show fairly large linewidths 
comparable to that in La 2 - a; Sr a ;Cu04. As it will be 
demonstrated below, the NQR data imply the same de- 
gree of the hole density inhomogeneity in La2- a: Sr :I ;Cu04 
and HgBa2Cu04+5. We will also show that despite hav- 
ing similar amplitudes of the charge inhomogeneity, the 
spatial profiles of the density distribution are very dif- 
ferent in these two cases: it is much smoother in the 
oxygen doped HgBa2Cu04+a. Consequently, the for- 
ward scattering is predominant and the mean free path 
in HgBa2Cu04+5 is substantially longer than that in 
La2- 2; Sr 2 ;Cu04. Our simulation also indicates that the 
screening of the Coulomb repulsion between the charge 
carriers in HgBa 2 Cu04 + 5 is stronger. We argue that 
these two reasons might explain the difference in super- 
conducting critical temperatures between these two com- 
pounds. 

The nanoscale charge density inhomogeneity in under- 
doped cuprates is an indirect way to distinguish the large 
Fermi surface of a normal Fermi liquid and the small 
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Fermi surface of a doped Mott insulator. Calculations in 
Ref. [h] for La2_ a: Sr a; Cu04 were based on the small hole 
pocket Fermi surface, which implies a small number of 
mobile charge carries. The small number of charge car- 
ries results in poor screenin g 10 i 13 and hence in the strong 
charge density inhomogeneity consistent with NQR data. 
On the other hand, the large Fermi surface implies the 
large number of mobile charge carriers and a very effec- 
tive Coulomb screening. In this case calculations give a 
very moderate charge inhomogeneit y 14 ' 15 , which is too 
weak to explain the observed NQR linewidths. 

Another important point is that the charge inhomo- 
geneity in the surface layer is always stronger than that 
in the bulk layer. This is because the surface ionic di- 
electric constant is about two times smaller than that in 
the bulki£. Hence the screening at the surface is weaker 
and this results in the stronger charge inhomogeneity 

In this paper, we present Hartree-Fock calculations in 
undcrdoped cupratcs with realistic parameters to simu- 
late the bulk charge distributions in HgBa2Cu04+5 and 
I^-zSr^CuOi, and the surface charge distribution in 
Bi2Si"2CaCu208+,5. Our calculation covers from under- 
doped to optimally doped regime. We fine tune theoreti- 
cal parameters to reproduce experimentally known NQR 
linewidths^i This enables us to perform a very accu- 
rate comparison of HgBa2Cu04+a and La2- x Sr :E Cu04. 
We find the charge inhomogeneity of a similar scale in 
both cases. However, the landscapes of spatial modu- 
lations and the hole mean free paths are substantially 
different. This is due to different positions of the dopant 
oxygen and the Sr-ion relative to the Cu02 plane. In ad- 
dition, the dielectric constants in these two compounds 
arc different. 

Motivated by STM dataSri in Bi 2 Sr2CaCu 2 8+< 5, we 
calculate also the hole density distribution in the surface 
Cu02 layer of this compound. The results obtained show 
a large charge density inhomogeneity comparable to that 
in the bulk of HgBa2Cu04_)- ( 5. Naturally, the local charge 
density is highly correlated with the interstitial oxygen 
positions as has been noticed previously^. 

Structure of the paper is the following: in Sec. II we 
formulate the effective model for an isolated Cu02 layer. 
Because of the long-range nature of the Coulomb inter- 
action, however, the isolated layer approximation is not 
sufficient and one has to take into account other layers. 
The impact of the other layers depends on the way of 
doping and on structural details: whether this is a sin- 
gle layer or a double layer compound. The single layer 
La2-xSr x Cu04 and HgBa2Cu04+a arc considered in Sec- 
tion III. In Sec. IV we simulate the charge distribution 
on the surface of the double layer Bi 2 Sr2CaCu208+«5. We 
also calculate the correlation function between the local 
charge density and oxygen dopant positions, and com- 
pare the result with the local DOS correlation function 
measured by STM. Our conclusions are presented in Sec 
V. 



II. THE EFFECTIVE MODEL FOR A Cu0 2 
PLANE 

We adopt the effective model formulated in Ref. [l^ 
based on the picture of a lightly doped Mott insulator. 
Throughout the paper we denote the average hole con- 
centration by p and assume p< 1, The central point is 
that the number of charge carriers is p instead of I — p as 
one would expect for a normal Fermi liquid. We consider 
first an "isolated" Cu02 layer. The "isolated" means 
that we disregard screening effects of the other layers. 

To construct a model relevant to cuprates, we first no- 
tice that there are the following distinct length scales: 
(i) The scale of the order of 1-2 lattice spacing where the 
doped holes are dressed by multiple virtual magnons. (ii) 
The scale about average separation between Coulomb de- 
fects or average separation between holes ~ 1/ ■ s Jp. This 
is the scale of the Coulomb screening. (Hi) The scale 
r ^> 1 / \fv~- The Coulomb gap may develop at this scale 
due to Anderson localization effects. 

Regarding the first point, we do not treat the strong 
correlations explicitly, but instead adopt the effective 
hole dispersion after quantum fluctuations at short dis- 
tances are included. To stress this point, we frequently 
use the term "holon" instead of "hole". It is known 
that the holon dispersion has minima at momenta ko = 
(±7r/2, ±7r/2), and is approximately isotropic around 
these points^. The bandwidth of the holon is about 
2J, where J « 130 mcV is the superexchange in the 
t-J model, although we do not directly employ the t-J 
formalism. Hereafter we set the following energy and 
distance units 

J = 130 meV 1 , 
ao = O.SSnm — > 1 , (1) 

where ao is the lattice spacing of the Cu02 plane. To 
imitate the holon dispersion we consider spinlcss fcrmions 
on a 2D square lattice. The Hamiltonian reads as follows: 

H t = Y,t"4c^ (2) 

m 

where cj is the holon creation operator at site i, and t" 
denotes the next-next-nearest-neighbor hopping on the 
square lattice. The Hamiltonian @ yields the following 
dispersion 

£ k = 2t"(cos 2k x + cos 2k y ) . (3) 

The dispersion is isotropic around minima at ko = 
(±7r/2, ±7r/2) as shown in Fig.Q] We choose t" = 0.25 J 
to reproduce the realistic holon band width, about 2 J, as 
obtained from numerical simulations of the t — J models 
An additional argument supporting this value is as fol- 
lows. Near its minimum the dispersion ([2]) can be ex- 
panded as =— > const + 4t"|k — ko| 2 , so the holon 
effective mass is equal to m* = Ti 2 / (Sa^t"). The value 
t" = 0.25J results in the effective mass m* ~ 2m e 
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FIG. 1: Dispersion minima of the spinless fermion generated 
by Hamiltonian @. 



which is close to the effective mass measured in mag- 
netic quantum oscillation experiment s] 1 8 1 19 , The realistic 
holon band width and the realistic effective mass jus- 
tify our choice of t" = 0.25 J. We note also that in the 
original t-J model formalism, there are four holon half- 
pockets inside magnetic Brillouin zone, and each pocket 
has two pscudospinsii; in the present model, we con- 
sider four full pockets inside the full Brillouin zone with 
spinless fermions, hence the number of charge degrees of 
freedom is exactly the same3 

Even though the holes are heavily dressed by magnetic 
fluctuations, their charge is conserved and hence they 
interact with Coulomb defects via the ordinary Coulomb 
potential: 



h-O 



Uu 



ZR 



+ A 2 



(4) 



Here is position of the holon and Ri is the in-plane pro- 
jection of the Coulomb defect (Sr-ion or dopant oxygen) 
position. The distance from the plane to the defect is A, 
and azB, ~ 0.8 is the size of the Zhang-Rice singlet. (We 
recall that the energy and distances are given in units of 
J and ao, correspondingly). The dimensionless "charge" 
Q ~ 0.5 depends on the compound, and we discuss its 
precise values later. 

Holon-holon Coulomb interaction is of a similar form 



Hi, 



Uijc\ciCjCj , 



many-body problem with the Hamiltonian (j6]) using the 
standard Hartree-Fock method. In other words we use 
the Hartree-Fock decomposition of the Coulomb interac- 
tion between holons: 



Hi. 



E^ 



'Ac, 



E 



U. l3 {c\ Cj ) 



c) Cl 



(7) 



This can be done for zero as well as for finite tempera- 
tures, as it is described in Ref. [IFj. 

The above formulation would solve the Coulomb prob- 
lem for the "isolated" C11O2 layer. However, the layer is 
always embedded in a multilayer structure, and because 
of the long-range nature of the Coulomb interaction, we 
have to take into account other layers. Their influence 
depends on the lattice structure. In the next section, we 
consider two different families of single layer cuprates. 



III. SINGLE LAYER HgBa 2 Cu0 4+ « AND 
La^Sr^CuCU COMPOUNDS. CHARGE 
DENSITY DISTRIBUTION, NQR LINE SHAPE, 
DENSITY OF STATES 



We treat a particular CuC>2 plane using the Hartree- 
Fock (HF) method. The role of other Cu02 planes is 
to provide screening of the Coulomb interaction in the 
Hartree-Fock analysis. Since CuC>2 planes have a very 
high longitudinal polarizability we consider the "other 
planes" as purely metallic. It has been shown in Ref. [h] 
that this "metallic approximation" is valid at p > 1 — 2% 
when the polarizability is sufficiently high. Within this 
approximation the HF plane in a single layer compound 
is sandwiched between two "metallic" sheets, as demon- 
strated in Fig. [5]for La2- 2; Sr 2 ;Cu04. 
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(5) 



"hh 



where a\ h 



^ a ZR 



1 stands for the combined size of 
two Zhang-Rice singlets. 
The total Hamiltonian 



H = H t + H, 



h-O 



Hi, 



(6) 



describes the in-plane Coulomb problem. Since the 
Coulomb interaction is not very strong, we solve the 
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FIG. 2: (Color online) The Hartree-Fock model of 
La2-a;Sr x Cu04 with four (m = 1,2,3,4) layers of Coulomb 
defects, i.e., doped Sr-ions, and with two "metallic" sheets. 
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A. La2- a; Sr a ;Cu04 

Coulomb defects in L^-xSr^CuC^ are created by Sr 
substitution for La ions. Therefore each defect donates 
one hole 



p = x 



(8) 



Given the periodic structure of CuC>2 planes along 
c— axis, there are two layers of Sr defects between two 
neighboring CuC>2 planes, as shown in Fig. [2] We mark 
the defect layers by the index m = {1,2,3,4}, so the 
HF CuC>2 plane is under the influence of 4 layers of Sr 
defects. Concentration of defects in each defect layer 
is x/2. The planar positions of defects are assumed to 
be random, and defects are assumed to sit above the 
center of the Cu plaquette. In addition, in each defect 
layer we impose a condition that defects never sit next 
to each other, i.e., the distance between defects is always 
larger than s/2. The distance between CuC>2 layers is 
d = 1.75, and the geometric distances to Coulomb defects 
are Ai = A4 = 1.15, A2 = A3 = 0.6. After accounting for 
the screening by "metallic" planes, the holon interaction 
with a defect (U]) is replaced by 



-Q 



E 



1 



y/\Ri -r,\ 2 + a dm \ 2 



(9) 



\ v /|R,-r,| 2 + (2nd+A, 



E 



(-1)' 



[ A /|R,-r i | 2 + (2nd- A m ) 2 / 

The summation over n reflects the image method for 
metallic screening^ With values of the distances A given 
above, the parameter a<j, defined in Eq. ((4]), takes the fol- 
lowing values: 0^/1 = 0^4 = 1.4, and a$x = fld3 = 1. The 
effective charge Q is determined by the dielectric con- 
stant of the lattice ei. In our previous workifi we used 
the value ei = 40 that approximately corresponds to the 
undoped compound.— This gave a reasonable fit of the 
NQR lines, but still the widths were about 20-30% larger 
compared to experiment. In the present paper, we fine 
tune the NQR widths by using q as a fitting parameter. 
Our best fit yields 



LSCO: ei = 30(1 + 6.25p) 



(10) 



Thus the value of e; depends on doping, e; = 30 for the 
undoped compound^, and e/ = 60 for the optimally 
doped compound (p = 0.16). The coefficient 6.25 has 
been obtained by fitting the NQR widths, see below. It 
is natural to have some doping dependence of the lat- 
tice dielectric constant since the lattice dynamics may 
change locally (get softer) with La — > Sr substitution. In 
our dimcnsionlcss units, the effective charge Q is 



Q 



eiaoJ 



30/e z 



(11) 



According to the same logic the Coulomb interaction 
between holons ([5]) is replaced by 



U, 



Q 



1 



r i\ + a HF 

2(-l)" 



(2ndy 



(12) 



The HF Hamiltonian is diagonalized in a 36 x 36 clus- 
ter. The resulting charge densities for particular realiza- 
tions of Coulomb defects at different doping levels are 
shown in Fig. [31 To calculate the NQR spectrum of in- 



p=0.02 




FIG. 3: The charge density of mobile holes in La2-xSr ;r CuO'4 
calculated at zero temperature and different doping levels. 

plane 63 Cu, we first calculate the hole density distribu- 
tion averaged over 20 random impurity configurations, 
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and then convert the hole density at a site i to the NQR 
frequency using the following scaling^ 



33 + 19m (MHz) 



(13) 



which implies that the NQR spectrum is directly related 
to the charge density distribution. The resulting NQR 
spectra for zero temperature and for T = 600 K are 
presented in Fig. @] The theoretical NQR spectra at 




FIG. 4: (Color online) Theoretical NQR spectra for in-plane 
6:! Cu in La2-xSr x Cu04 at T = and T = 600 K and for 
different doping levels p. 

T = 600 K agree very well with experimental data 6 . 
Although we are not aware of NQR data at low tem- 
peratures, the theoretical spectra at T = are shown in 
order to demonstrate how the charge density distribution 
evolves with temperature. It is instructive to present the 
widths of the charge density distribution. Our simulation 
yields 



LSCO at p = 0.16 



LSCO atp = 0.08 



r P (T 
r 



P (T-- 



0) = 0.12 
600*0 = 0.09, 



T P (T= 0) = 0.07 
T P (T= 600*0 = 0.05. 



(14) 



In our previous workjiS we performed similar calcu- 
lations for La2- K Sr :E Cu04 and demonstrated that the 
Hartree-Fock calculation reproduces reasonably well the 
evolution of NQR lincshapes with doping. However, the 
calculated NQR lincwidths were about 20-30% larger 
compared to the experimental values. The goal of the 



present calculation is to fit the NQR lincwidths quan- 
titatively, using the dielectric constant (|10p as a fit- 
ting parameter. Based on this fit we accurately quan- 
tify the charge inhomogeneity and the hole mean free 
path as discussed below. We also note that the previ- 
ous calculation^ has found an additional high frequency 
hump/shoulder in the NQR line. The hump was due to 
the strong Coulomb binding of holes to the accidentally 
clustered Coulomb defects. The binding is very sensitive 
to the strength of the Coulomb attraction. This binding 
does not show up in the present calculation because of 
the larger value of the dielectric constant used, and also 
due to the fact that we impose here the "nonadjacent" 
condition for Coulomb defects [see the paragraph before 
Eq. ©]. As suggested in Ref. [a, the experimentally ob- 
served high frequency hump (the NQR " B-line" ) is most 
probably due to the direct action of Sr ion Coulomb po- 
tential on the Cu nucleus, which is beyond the scope of 
the present model. 



B. HgBa 2 Cu0 4+a 



Interstitial oxygen ions in HgE^CuO^a are located 
right between two neighboring Cu02 layers, on top 
of four adjacent Cu sites^i The minimal model for 
HgBa2Cu04+a is shown in Fig. [5j There is one layer of 
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FIG. 5: (Color online) The Hartree-Fock model of 
HgBa2Cu04+i with two layers of Coulomb defects (intersti- 
tial oxygen ions) and two "metallic" sheets. 



oxygen defects above the Cu02 plane and another defect 
layer below it. Concentration of defects in each defect 
layer is 8. We assume that each interstitial oxygen do- 
nates two holes into the Cu02 plane, i.e., 



p = 26 



(15) 
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According to this picture, the holon interaction with a 
defect, Eq. ^ is replaced by 



-2Q 



E 



(16) 



(-1)" 



^ V|R/-r,| 2 + (2?irf + d/2) 2 



(-l) r 



+ Ji A /|R i -r i p + (2nd-d/2)2 y / " 

The interlayer distance in HgBa2Cu04+,5 is rf = 2.5, 
and each interstitial oxygen carries "charge" — 2Q. Tak- 
ing into account size of the Zhang- Rice singlet, we have 
ad = \J a 2 ZR + (d/2) 2 sa 1.5. The hole-hole interaction 
has exactly the same form as that in La2— xSrxCuC^, 
Eq. (|12|). with the interlayer distance d = 2.5 and with 
the value of Q described below. In each defect layer we 
simulate positions of the defects randomly, and impose 
again the nonadjacent condition for the defects. 

Similar to the procedure in La2- x Sr T Cu04, we fine 
tune the NQR widths using the lattice dielectric con- 
stant ei as fitting parameter. We find that in order to 
fit experimental NQR spectra^iii the dielectric constant 
has to be taken as 



et = 30(1 + 25p) 



(17) 



Note that the doping dependence of e/ in this case is 
four times stronger than that in Eq. (JTUJ). The stronger 
dependence is quite natural since it is due to the shift 
of the interstitial oxygen position in an applied electric 
field. The shift is significant because binding of inter- 
stitial oxygens in the lattice is weak. In other words, a 
dopant oxygen brings in new local phonon modes which 
enhance the dielectric constant. The effective charge Q is 
determined by the same Eq. (jlip . with e; from Eq. (|17p . 

The HF Hamiltonian is diagonalized in a 36 x 36 clus- 
ter. The resulting charge densities for particular realiza- 
tions of Coulomb defects at different doping levels are 
shown in Fig. [51 A strong charge inhomogeneity in- 
duced by oxygen dopants is apparent. However, den- 
sity distribution profiles here are much smoother than in 
I^-zSr^CuC^ as one can easily see by comparing Figs. [6] 
and [3] at same doping levels. 

Now, we address the NQR spectrum of the in-plane 
63 Cu in HgBa2Cu04+5. We first calculate the hole den- 
sity distribution averaged over 20 random impurity con- 
figurations, and then convert the hole density at a site i 
to the NQR frequency by 



C + 30n, (MHz) 



(18) 



where C is a doping independent constant. Note that 
the coefficient 30 MHz in this formula is different from 
that in Eq. ([T3"|). We found the coefficient in Eq. ([IS]) by 
fitting the NQR-line centers using the experimental data 
of Ref. [ll] (this way is more reliable than using results 




FIG. 6: The charge density of mobile holes in HgBa2Cu04+,5 
at zero temperature and at different doping levels. 
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FIG. 7: (Color online) Theoretical NQR spectra for in-plane 
63 Cu in HgBa2Cu04+« at zero temperature and for different 
doping levels p. 



of theoretical calculations of electric field gradients^) 
Since NQR data for HgBa 2 Cu04 +( 5 are taken at suffi- 
ciently low temperatures, we perform calculations only 
at T = 0. Fig. [7] shows the theoretical NQR lines by 
assuming O = 0, since it is an irrelevant constant shift. 
The calculations agree well with the experimental data 
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of Refs. 0IH At optimal doping p = 0.16, the NQR 
linewidth in HgBa 2 Cu0 4+( 5 is T NQR w 2.6 MHz. The 
widths of the charge density distribution arc 

HBCO atp = 0.16 : T p (T = 0) = 0.09 , 
HBCO at p = 0.08 : T P (T = 0) = 0.08 . (19) 

Comparing with the corresponding values for 
La2_ a; Sr a ;Cu04 in Eq. (fT4|) . we see that at zero 
temperature the amplitudes of charge inhomogeneities 
in HgBa2Cu04-|-5 and in La2- x Sr x Cu04 are quite 
similar. At optimal doping p = 0.16, the amplitude of 
charge inhomogeneity in HgBa 2 Cu04 + a is only by 30% 
smaller than that in La 2 _ x Sr :r ,Cu0 4 . 



C. Density of states and mean free path 

A comparison between Eqs. (|14[) and ([TO)) suggests 
that the hole density inhomogeneity in HgBa2Cu0.4-|_ < 5 is 
pretty close to that in La2- a jSr a: Cu0 4 . This seems rather 
surprising, given that the superconducting critical tem- 
perature in HgBa2Cu04+5 is much higher than that in 
La2_ 2; Sr 2 ;Cu04. Indeed, a disorder reduces the d— wave 
superconducting critical temperature and, based on this 
argument, one can expect that La2- a; Sr a ;Cu0 4 is "more 
disordered" . However, one should quantify exactly what 
is the measure of the disorder. While the overall ampli- 
tude of the charge inhomogeneity is one possible factor, 
the hole mean free path is another and in fact more im- 
portant measure of the disorder. We already noticed that 
the spatial charge distribution profiles induced by inter- 
stitial oxygens and Sr-dopants are rather different. To 
quantify this difference in more detail, we calculate the 
quasiparticle DOS defined in the standard way 

p(e) = ^J2 s ^- E -^ ( 2 °) 

n 

where N is the total number of sites, and E n is the n-th 
eigenenergy. Plots of the DOS in HgBa 2 Cu04-|- < $ and in 
La2- 2; Sr :Z ;Cu04 at doping p = 0.16 are shown in Figs. [5] 
and [9l together with the charge profiles for particular 
disorder realizations. The DOS has been obtained af- 
ter averaging over 40 different disorder realizations. We 
have checked that the DOS calculated with 20 realiza- 
tions of random dopant positions is practically the same 
as the DOS calculated with 40 realizations. This means 
that 40 realizations is sufficient to disregard the statisti- 
cal noise. The calculated DOS exhibits pronounced os- 
cillations. These oscillations arc a byproduct of the fi- 
nite size of the cluster. Maxima of the DOS correspond 
to degenerate states with dispersion ([3]) on the 36x36 
torus. The oscillations must certainly disappear in the 
thermodynamic limit. However, oscillations have a phys- 
ical meaning: they indicate that the quantum states are 
quite extended, with the hole mean free path I compa- 
rable with the size of the cluster used. To estimate the 




FIG. 8: (Color online) (a) Density of states and (b) charge 
density profile in HgBa2Cu04+,5 at optimal doping p = 0.16. 
Shown in the inset of panel (a) is the DOS of the homogeneous 
model with effective scattering rate r] — 0.014. 




FIG. 9: (Color online) (a) DOS and (b) charge density profile 
in La2- a: Sr a; Cu04 at p = 0.16. The inset shows the DOS of 
the homogeneous model with scattering rate r\ — 0.021. 
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mean free path more accurately we use the following pro- 
cedure. We consider the DOS of an ideal homogeneous 
system described by the Hamiltonian ([2]). It consists of 6- 
functions whose positions are fixed by periodic boundary 
conditions. The weight of each (5-function can be eas- 
ily calculated for the 36x36 torus. Now, we artificially 
broaden each (5-function, 



S(e - e„) 



1 



V 



7T (e - e„) 2 + rj 1 



(21) 



to simulate a disorder scattering. We find that the DOS 
of this model is very sensitive to r\. Then, we adjust 
the broadening rj to reproduce (roughly) the amplitudes 
of the DOS oscillations obtained in our actual calcula- 
tions, as shown in Figs. [Ha) andO[a). This gives the 
following values of the effective broadening: r\ = 0.014 
for HgBa 2 CuO4.08, and rj = 0.021 for Lai.84Sr .i6CuO4. 
The mean free path of a hole at the Fermi surface can be 
then estimated as: 



I = VpT 



vf 
2 V ' 



(22) 



where r = l/2rj is the collision time. At small doping, 
p -C 1, the dispersion ([3]) results in the Fermi velocity 
Vf ~ 8t"aoy/irp. Hence, at p = 0.16, the Fermi velocity 
is vf ~ IAJoq. Together with the above values of 77, this 
results in the following estimates for the hole mean free 
paths: 



HgBa 2 CuO im : 
Lai. 84 SV . i6 Cu0 4 



I 



50ao : 
34a 



19nm 
linm 



(23) 



We thus find that La2- a: Sr :l: Cu04 is indeed more disor- 
dered, in the sense that it has a shorter mean free path. 
The different mean free paths are due to different dis- 
tances from the Cu02 plane to the Coulomb defect (in- 
terstitial oxygen versus Sr-dopant), and due to different 
distances between Coulomb defects. In HgBa2Cu04+,5, 
both distances are larger, therefore the Coulomb poten- 
tial is smoother and hence less contributes to the large 
angle scattering. This difference is clearly reflected in the 
charge density profiles shown in Fig.HJb) and^b). Over- 
all amplitudes of the charge inhomogeneity are pretty 
close, but the charge distribution in La2- a: Sr a: Cu04 is 
much more spiky. This gives rise to the strong scattering 
with large momentum transfer which is detrimental for 
d— wave superconductivity. 

Considering the superconducting correlation length 
£ ~ 2 nm, one finds that Z/£ ~ 6 in La2- a: Sr :l: Cu04 
and Z/£ ~ 10 in HgBa 2 Cu04+,5. One should there- 
fore expect that the disorder suppression of supercon- 
ductivity in HgBa 2 Cu0 4 +,5 is indeed weaker. In addi- 
tion, according to our fits of the NQR data, the effective 
ionic dielectric constant ei in HgBa 2 Cu0 4+ 5, Eq. (fl~7|) . 
is larger than that in La2_ 2; Sr 2 ;Cu04, Eq. (jTU|) . At the 
optimal doping p = 0.16, for example, ei(LSCO) = 60 
and ei(HBCO) = 150. The larger e; in HgBa 2 Cu0 4+ «5 



implies better screening of the Coulomb repulsion be- 
tween holes and hence a smaller Coulomb pseudopoten- 
tial. This may further enhance T c of HgBa2Cu04+,5 com- 
pared to that of La2- 2: Sr ;r Cu04. 



IV. NANOSCALE HOLE DENSITY 
INHOMOGENEITY IN THE SURFACE LAYER 

OF BiaSraCaCuaOg+i 

Recent STM experiments^— have revealed large vari- 
ations of pairing gaps in Bi2Sr2CaCu20s+5 which are 
spatially correlated with the dopant oxygen density. Mo- 
tivated by this observation, we now consider the charge 
distribution in the surface Cu02 layer of this compound. 
Our purpose is to estimate the magnitude of the charge 
inhomogeneity induced by interstitial oxygen dopants, 
and see how the hole density profile is correlated with 
the position of these oxygens. In the present work, wc 
do not calculate the local DOS and hence cannot address 
STM directly. Instead, we calculate the charge distribu- 
tion similar to the previous sections, and show that the 
effect of interstitial oxygen on charge inhomogeneity is 
very significant. 

As we have seen above, charge distribution profiles 
depend on the lattice dielectric constant whose val- 
ues can be reliably obtained by fitting the NQR data. 
Due to lack of systematic NQR data, we cannot de- 
termine the dielectric constant of Bi2Sr2CaCu 2 08+,5 in 
this way. Instead, we assume that the doping de- 
pendence of the lattice dielectric constant in the bulk 
of Bi2Sr2CaCu208+5 is described by the same formula 
Eq. (fT7|) as in HgBa 2 Cu04 + 5. This is because both com- 
pounds arc doped by interstitial oxygens which have sim- 
ilar influence on lattice dynamics and dielectric screen- 
ing. This adoption certainly ignores lattice structure de- 
tails and should be further refined when NQR data is 
available. Nevertheless, we found that our results for 
Bi2Sr2CaCu 2 08+<5 are fairly robust within a sensible vari- 
ations of the dielectric constant, and thus they should 
give a qualitative description of the surface charge inho- 
mogeneity. 



surface 



BiO 



HF 



metal 



-2Q 



d 



FIG. 10: (Color online) The Hartree-Fock model for 
BiaSraCaCuaOg+a surface. Blue dots indicate the interstitial 
oxygen dopants. 

Bi 2 Sr2CaCu208+<5 is a double layer compound and a 
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schematic picture of the two surface CuC>2 layers is shown 
in Fig. [TU1 Interstitial oxygens are assumed to be lo- 
cated at A w 0.37 nm above the CuC>2 layer, which is 
shown to be the most energetically favorable position^ 
We assume that the surface concentration of Coulomb 
defects (interstitial oxygens) is 5. They are randomly 
distributed and, similarly to the previous considerations, 
we impose again the condition that the defects never sit 
next to each other, i.e., the distance between the defects 
is always larger than y2 (in units of do). The separa- 
tion between the CuC-2 layers in the double layer struc- 
ture is d = 0.33 miK 0.86. We treat the top layer by 
Hartree-Fock approximation, and the underneath layer 
as a "metallic" sheet that provides screening. There are 
also Coulomb defects underneath of the screening layer, 
but they are well screened by the metallic sheet and hence 
do not influence the Hartree-Fock procedure. The hole- 
defect and the hole-hole interaction in the Hartree-Fock 
top layer are of the following form: 



p=0.08 



-2Q 




Ui 



v /|Ri-r l | 2 + (2rf + A) 2 / 
1 




(24) 



Here = -y/A 2 



1.28, and cihf 



ZR 



1. 



Note that there is only one image charge per physical 
charge because there is only one screening layer. We as- 
sume that the lattice contribution to the dielectric con- 
stant in Bi2Si'2CaCu208+5 has the form Eq. ([T7|) . On the 
surface, the dielectric constant is expected to be reduced 
to half of its value in the bulk,— therefore 



^BSCCO 



15(1 + 25p) 



(25) 



and the value of effective "charge" Q follows from 
Eq. (UU). 

Plots of the in-plane hole density for 8 = 0.04 (p = 
0.08) and 5 = 0.08 (p = 0.16) are presented in Fig. HQ 
together with corresponding maps for these particular re- 
alizations of Coulomb defects. The plots demonstrate a 
large spatial variations in hole density. For the S = 0.04 
case, the range of hole density modulation is about 
0.03 < n, < 0.15, which is a very significant fluctua- 
tion, having in mind that the average hole density is 
(jn) = p = 0.08. The same strong inhomogeneity is 
seen in the 6 = 0.08 case, where the local density varies 
in the range roughly 0.10 < n.; < 0.23, while average 
density is = p = 0.16. The density distribution 
curves for different doping levels are shown in Fig. 1121 
We stress that the precise density profiles depend on the 
dielectric constant which, due to lack of NQR data for 




(c) 



(f) / "» 



FIG. 11: (Color online) (a,d) The zero temperature hole den- 
sity plots and (b,e) the hole density maps for the surface layer 
of Bi2Sr2CaCu2 0-8+6 for the average hole density p = 0.08 
(left panels) and 0.16 (right panels). The lower panels (c) 
and (f) show positions of the oxygen dopants for the corre- 
sponding realizations of disorder. It is evident that the dopant 
oxygens locally increase the hole density. 



Bi2Sr2CaCu2 Os+a, is taken here in an ad hoc way. On 
a qualitative level, however, the density inhomogeneity 
is rather stable with respect to the value of the dielec- 
tric constant. For example, the widths of the surface 
density distributions plotted in Fig. Q2] are only slightly 
larger than those obtained in the case of HgBa2Cu04+,5, 
see Eq. (TT9|l . This is in spite of the fact that the re- 
sults for HgBa2Cu04+a have been obtained at the twice 
larger value of the dielectric constant. Therefore, we be- 
lieve that our conclusion about the strong surface charge 
density inhomogeneity is very reliable. 

Naturally, the charge density of mobile holes is higher 
around areas with higher interstitial oxygen concentra- 
tion. Fig. [TT] clearly shows this correlation. To quantify 
this, we define a correlation function between local hole 
density and dopant oxygen positions, which is analogous 
to the correlator introduced in Ref. for analysis of the 
spatial variations of the local DOS. On a discrete lattice, 
the hole density rij is defined on sites of the square lattice 
representing the Cu02 plane. The function fi indicates 
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FIG. 12: (Color online) The zero temperature hole density 
distribution in the surface layer of Bi2Sr2CaCu208+,5 at var- 
ious doping levels p. Each curve is averaged over 20 realiza- 
tions of disorder. 



the location of interstitial oxygens: 



fi = 



1 if i € dopant oxygen 
elsewhere , 



(26) 



where i runs through points at the center of the pla- 
quettes. The density-oxygen correlation function is then 
defined as 



, s 1 \ - [Si - 7] [rtj+R - n] 
C »'W-n\ ^A n 

with proper normalizations 

i 



(27) 



(28) 



where / = j?J2iSi = s and n = jf J2j n i = V- Fig. EH 
shows the correlation function C„/(R) averaged over 20 
disorder realizations for each doping. There is a clear 
positive correlation due to the Coulomb attraction to the 
oxygen defects. The value of the correlator at R — > is 
C ~ 0.3 — 0.4 and the scale at which it goes to zero 
is about 10-15 A. Interestingly, the correlator between 
the gap in local DOS and the interstitial oxygen position 
measured in STM shows the same positive correlation 
with very similar scales^. Further investigation is nec- 
essary to clarify if there is a connection between these 
two correlators, and to understand the physical reasons 
behind this apparent correspondence. 



In this paper, we study the spatial distribution of 
doped holes in cuprates, focusing particularly on a 
comparison between two different physical situations: 
doping by a cationic substitution Sr for La as in 
La2- 2; Sr 2 ;Cu04, and doping by interstitial oxygen ions 
as in HgBa2Cu04_|-5. The main results are summarized 
as follows. 

The hole density inhomogeneity in HgBa2Cu04 + a is 
nearly as strong as in La2- K Sr 2 .Cu04. For example, at 
the optimal doping p = 0.16 the width of the charge den- 
sity distribution is about T p = 0.09, which is close to the 
corresponding number T p = 0.12 in L^-ajSr^CuC^. This 
conclusion is well supported by the comparison of our cal- 
culations with the existing NQR data in HgBa2Cu04+5 
and La2- x Sr x Cu04. In spite of the close overall ampli- 
tudes, the landscape of charge inhomogeneity in these 
two compounds are very different. In oxygen doped 
HgBa2Cu04+5, the disorder potential profiles are much 
smoother than that in Sr-dopcd L^-zSr^CuC^. Corre- 
spondingly, the hole mean free path in HgBa2Cu04+5 
is larger. In other words, disorder induced scatter- 
ing processes with a large momentum transfer (which 
are destructive for d— wave pairing) are less pronounced 
in oxygen doped HgBa2Cu04 +< 5 compared to the case 
of La2_ x Sr :E Cu04. In addition, the screening of the 
Coulomb repulsion between holes in HgBa2Cu04+,5 is 
about twice stronger than that in La2- :E Sr. E Cu04. In 
our opinion, these two reasons might explain the much 
higher superconducting critical temperature of oxygen 
doped HgBa2Cu04 +< 5. 

We found that the charge density nanoscale inhomo- 
geneity in the surface Cu02 layer of Bi 2 Sr 2 CaCu208+5 
(the layer available for STM) is of the same magnitude as 
that in the bulk of HgBa2Cu04+5. As expected on phys- 
ical grounds, the hole density positively correlates with 
the positions of interstitial dopant oxygen. Remarkably, 
the correlation function obtained here resembles the pos- 
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itive correlation between the local gap and dopant oxy- 
gens seen in the STM data. The reason for this appar- 
ent coincidence and implications of the charge inhomo- 
geneity for the spatial variations of the pairing gaps in 
Bi2Sr2CaCu208+<5 have to be clarified in future studies. 
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